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Abstract. We consider a two-dimensional crystalline monolayer of carbon atoms with a single 
non-magnetic impurity. Using the Weyl Hamiltonian to describe electronic energy spectrum near 
the Dirac points, we calculate the wave function and energy of impurity states, as well as the 
induced magnetic moment associated with polarization of the electron system. We present a 
phase diagram of the localized magnetic state as a function of the chemical potential /i and 
coupling constant gc- 



1. Introduction 

Graphene has been attracting a lot of attention recently due to its unusual electronic structure 
and, correspondigly, due to its unique transport properties, which differ substantially from those 
of traditional metals and semiconductors [H El |3] . The key characteristic features of graphene are 
the exact two-dimensionality of its crystal structure and massless relativistic energy spectrum 
near the Fermi surface. The Fermi surface transforms into isolated points in the clean limit of 
undoped graphene, and the electronic spectrum in the vicinity of these points can be described 
by the relativistic Dirac model. It is now commonly believed that these unusual properties 
of graphene make it an excellent candidate for various device applications in nanoelectronics 
and/or spintronics. [4j. 

Although there is great interest in magnetic properties of graphene, these are still not well 
understood. The main problem is that the usual doping with magnetic impurities like Mn or Fe 
can be ineffective because the magnetic atoms change dramatically the structure of electronic 
states due to bonding with carbon atoms. However, the possibility of doping with magnetic 
impurities, which can lead to ferromagnetism in graphene, has been discussed in several recent 
publications |Sl El E] • The main interest now is related to the possibility of magnetism induced 
by the presence of non- magnetic defects jH El [TU] . 

In this work we consider formation of a local magnetic moment in graphene, induced by a non- 
magnetic impurity. Assuming exchange coupling between the localized spin and electronic states 
we calculate magnetic polarization of the system and analyze conditions for the appearance of 
a local magnetic moment. 



2. Model of graphene with a single impurity 

The energy spectrum of graphene is described by the fohowing tight- binding Hamiltonian [H] , 



'Htb = X! ^AC^jCAi + ^ SBC^Bn^Bn + t {c 



(1) 



i n <i,n> 



where the creation and annihilation operators correspond to electrons localized in atomic sites 
of two nonequivalent sublattices A and B, with the hopping between the nearest neighbors 
of different sublattices described by the parameter t. The indices i and n refer to sites in the 
sublattices A and B, respectively. The onsite energy in graphene is the same for both sublattices, 
£A = £b- 

By transforming Hamiltonian (1) (with ea = £b) to the momentum representation, one 
reproduces the energy spectrum of graphene, which has two nonequivalent Dirac points in the 
k space, with a linear energy spectrum near these points. The Dirac points are known to be 
the most interesting features of the electronic structure, because the Fermi level of graphene is 
usually located at these points or in their close vicinity. 

In our calculations we restrict ourselves to the energy spectrum near the Dirac points. The 
corresponding Hamiltonian is known in quantum field theory as the 2D Weyl Hamiltonian of 
massless relativistic particles. 



where v is the velocity of electrons, which in graphene is f ~ 10^ m/s, r are the Pauli matrices 
acting in the sublattice space (matrices of pseudospin) , and we put the energy onset in the Dirac 
points, Ea = £b = ^■ 



Figure 1. Schematic view of an impurity C, which substitutes a carbon atom in the sublattice 
A of graphene. 

In the following we assume that one carbon atom in the sublattice A is replaced by an 
impurity atom C (located at r = 0), as shown schematically in Fig. 1. This replacement can 
be formally included in the Hamiltonian (1) as a term with onsite energy ec at the site i = 0, 
so that the perturbation is Vq = [ec — £A)c^AQ^Aa- After reducing to the Weyl Hamiltonian, the 
corresponding perturbation term in the polar coordinates (p, (p) 







(3) 



3. Localized electronic states 

Our task now is to find impurity states described by the Hamiltonian Ti = Ti.Q + 'Hint, with Ho 
and TCint defined by Eqs. (2) and (3), respectively. The corresponding Schrodinger equation, 
{7i — 9?) = 0, gives us the energy and eigenfunction of the localized electronic states. 

Solving the Schrodinger equation for p 7^ we find a general solution in the form of a 
superposition of the Bcsscl functions Jm(0 ^^^^ where ^ = pe/fiv) and m is the magnetic 

quantum number. It turns out that only the following pseudospinor form of the solution gives 
us the localized states near the Dirac point: 



where A'' is the normalization factor. 

Energy of the localized states can be found by integrating the Schrodinger equation with 
the wavefunction (4) over a small area near the point p = 0. Within the lattice model, this 
corresponds to integration over p < uq, where oq is the lattice parameter. This way we account 
for the boundary condition at /> = with the 5- like perturbation localized at this point. The 
ground state corresponds to the choice m = Oorm=— 1. Both states give us the same energy 
of the localized state. 



As we see from the solution (5), the energy of the localized state eo ^ as the perturbation | Vo| 
increases. For the parameters typical of graphene, the sign of logarithm in (5) is negative, and 
therefore the sign of the energy state £0 is opposite to that of perturbation potential. In other 
words, the attractive (negative) potential, Vq < 0, gives us a state with the energy eo > and 
vice versa. This is just the opposite behavior as compared to the impurity state in traditional 
semiconductors - if we assume that graphene is similar to the gapless semiconductors. An 
important point is the above mentioned degeneracy of localized states, which is related to the 
pseudospinor wave function in graphene. 

We have also checked the result (5) for the energy of impurity state by calculating the 
exact T-matrix of scattering from the impurity described by the matrix perturbation (3). Both 
calculation methods give exactly the same results. 

4. Induced magnetic moment 

The impurity state in graphene is occupied by an electron when the Fermi level is above the 
impurity level (5). Hence, in our consideration we assume that the location of the Fermi 
level is an independent parameter which is not related with the single impurity state under 
consideration. In reality, this assumption is well justified because position of the Fermi level in 
graphene is usually related with various defects or simply can be controlled by a gate voltage. 

The single electron localized at the impurity has the magnetic moment pB- In the following 
consideration we assume that the Coulomb interaction does not allow for two or more electrons 
to occupy the same impurity state, which means that the Hubbard energy is sufficiently large. 

Magnetic moment of the localized electron polarizes the electron system, inducing magnetic 
moment whose magnitude depends on the magnetic polarizability of the electron gas in graphene. 
Assuming exchange coupling of electrons in graphene with the localized electrons we calculate 
the induced magnetic moment. The Hamiltonian of such an exchange interaction is 




(4) 



271 h^v^ 27rh^v'^ 



(5) 



£0 = 



Vb In (ao |£o I /^ii') Vo In (27rao hv/\Vo\) 



(6) 



where Qc is the couphng constant, /iz(r) is the spatial distribution of magnetization associated 
with the wavefunction profile of the localized electron in Eq. (4), 

/.,(r)=5/iB|^(r)|2, (7) 

and (Tz is the Pauli matrix for electron spin. We take the quantization axis along the spin 
orientation of the localized electron. 

The induced magnetic moment M(r) of electrons in graphene can be calculated using 
the quantum field theoretical method, used for calculation of magnetic polarizability 
Accordingly, one can present it as a loop Feynman diagram, which gives 

M,(r) = ig^g^iB TV J d\' J^^zGo {e, r - r') a,Go (e, r' - r) /i,(r'), (8) 

where Go(e,r) is the Green function of electrons in graphene in the energy-coordinate 
representation. The calculation of the Green function using the Hamiltonian (2) gives 

where Hq^I are the Hankel functions \12l,- Using Eqs. (7)-(9) we find 

^ ' A-nfiv J |r-r'P' ^ ' 

where 

oo 

c = j die[Ki{i)+Ki{i)\, (11) 



and 

„-m+l 

Km{i) = -^Hg\tO (12) 

is the the modified Bessel function |12] . 

Using Eqs. (10)-(12) one can write the induced magnetic moment in the following form: 

snv Jo {P + P'){p-P'y ^ ^ ^ 

where £ (?) is the elliptic integral of the second kind and <j = 2-^ pp' /(p + p'). This expression 
was used for the numerical calculations, and the results are presented in Fig. 2. 
We have also calculated the total induced magnetic moment 

00 

MoiVo,gc) =2tt [ M,{p) pdp (14) 







The dependence of Mq on the strength of impurity potential |Vo| is found to be rather weak, 
whereas the main factor which may enhance Mq is the magnitude of the coupling parameter gc- 
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Figure 2. The density of induced magnetic moment as a function of the distance p for different 
couphng parameters. 

5. Magnetic coupling of the localized spin with induced magnetic moment 

The magnetic interaction of the locahzed spin with induced magnetic density Mz{p) leads to 
the renormahzation of the impurity energy level [IS] . This interaction can be written as 

oo 

Eint = -27rg,J M,{g)\<if{Q)\^ pdp. (15) 


We have calculated numerically the interaction (15) and found the corresponding shift of the 
energy level. Our results are presented in Fig. 3. 
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Figure 3. Energy level of the localized electron interacting with the induced magnetic moment 
for different values of the coupling constant gc- The energy in the units of hopping t. 




Figure 4. The phase diagram with the dashed area corresponding to the local magnetic moment 
induced by a nonmagnetic impurity. 

6. Phase diagram 

The obtained results can be presented in the form of a phase diagram, Fig. 4, which shows 
the regions of the parameters Vq and /U, for which the impurity level is occupied with electron 
inducing the local magnetization around impurity. It should be pointed out that the level can 
be occupied even if the nonrenormalized position of the level is above the Fermi energy |13j . 
This result is related to the effect of onsite magnetic correlations. 

7. Conclusion 

The results of our calculations show that the nonmagnetic impurity can create an induced 
localized magnetic state due to the local electronic correlations. It should be emphasized that 
this effect is not related to free electrons because the density of free electrons in graphene can 
be vanishingly small. However, the magnetic polarizability of the electron system is large and 
long-range even without free electrons, which is the main reason for the local ferromagnetism. 
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